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Abstract. Helium-rich subdwarfs are a rare subclass of hot subdwarf stars which con- 
stitute a small and inhomogeneous group showing varying degrees of helium enrichment. 
Only one star, LS IV ° 14 116 has been found to show multiperiodic luminosity variations. 
The variability of LS IV °14 116 has been explained as the consequence of nonradial 
g-mode oscillations, whose excitation is difficult to understand within the frame of the 
• Standard /c-mechanism driving pulsations in sdBV stars. In a recent study, we have pro- 

posed that the pulsations of LS IV °14 116 might be driven through the e- mechanism 
acting in unstable He-burning zones in the interior of the star, that appear before the qui- 
escent He-burning phase. One of the few accepted scenarios for the formation of He-rich 
subdwarfs is the merger of two He-core white dwarfs. As part of this project, we present 
a study of the e-mechanism in post-merger remnants, and discuss the results in the light 
of the pulsations exhibited by LS IV ° 14 116. 

1 Introduction 

About 5% of all hot subdwarf stars (sdB, sdO) show helium (He)-enriched surface abundances (He- 
sdB, He-sdO). While most normal hydrogen (H)-rich subdwarfs are supposed to be low mass core 
He-burning stars with atmospheres dominated by H due to the action of gravitational settling, the 
evolutionary status of the He-rich subclass is less clear. He-rich subdwarfs have been suggested to be 
the result of either double He white dwarf (HeWD) mergers [ 1] or late helium core flashes (a.k.a. hot- 
flashers[21). In fact, it is by no means clear whether intermediate and extreme He-rich sdB stars form a 
single class of objects or not Q. Regardless of the particular scenario proposed for their formation, it 
is accepted that some of these stars are still contracting towards the He-core burning phase (EHB), as 
otherwise gravitational settling of the remaining H would have turned the star into a H-rich sdB star. 

Two families of pulsators have been discovered within the H-rich sdB stars: the rapid pulsators 
(sdBVr; |4|) discovered by |5 1, and the slow pulsators (sdBVs; |4|) discovered by [6|. While sdBVr 
stars show short pulsation periods (~ 80 - 400 s) ascribed to radial modes and non-radial p-modes, 
sdBVs pulsations (with periods ~ 2500 - 7000 s) are associated to non-radial long period g-modes. 
Pulsations in both groups of variable stars have been explained by the action of the /c-mechanism due 
to the partial ionization of iron group elements in the outer layers, where these elements are enhanced 
by the action of radiative levitation |7| |8|. While many sdBV stars are known, only one He-sdB star, 
LS IV-14°116, has been found to be a pulsator Q ITO). In fact, LS IV-14°116 is a very intriguing 
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object: while its atmospheric parameters (Teff and g) place it well within the sdBVr instability region 
lfmifT2]| . it displays periods typical of the sdBVs family of pulsators fTO) lfT2]| . From a spectroscopic 
point of view, LS IV-14°116 is also an intriguing object, showing a mild He-enrichment (nee = 0.16) 
and very high excesses of zirconium, yttrium and strontium in its atmosphere [13 1. In this sense, both 
the He-enrichment and the heavy metal excesses might be related to the effects of ongoing diffusion 
before reaching diffusive equilibrium, thus placing the star in the pre-EHB phase. Up to now, the 
driving mechanism behind the long period pulsations of LS IV- 14° 116, as well as the absence of short 
period pulsations, remains a mystery and the star poses a challenge to the theory of stellar pulsations 
lfT4l . In this connection, the recent confirmation of both the multiperiodic variability and the T^ff - g 
values for LS IV- 14° 1 16 fl^ strongly increase the enigma. In a recent article |15 1 we have suggested 
that pulsations in LS IV- 14° 1 16 might be excited by the e-mechanism acting on unstable He-burning 
shells that take place before the begining of the EHB. In the e-mechanism, the driving is due to the 
strong dependence of nuclear burning rates on temperature. In the layers where nuclear reactions take 
place, thermal energy is gained at compression by the enhancement of nuclear energy liberation, while 
the opposite happens during the expansion phase 1 16|. In particular, in 1 15| we have shown that the 
e-mechanism is able to excite pulsations in hot-flasher stellar models. Although only one hot-flasher 
sequence was analyzed in ifTSl . it was shown that the e-mechanism is qualitatively able to reproduce 
the long period {P > 1000s) ^-modes at the values of \ogg and log red inferred for LS IV-14°116. 
As post- merger models also undergo off-centered He-shell flashes before settling on the EHB [1] and 
given that the location of the burning shell is particularly relevant for the range of excited periods fT5l . 
it seems natural to explore the effect of the e-mechanism in post He WD merger stellar models. 

Interestingly enough, no star has been identified so far to be excited by the epsilon mechanism. 
Aside from our previous suggestion! 15 1, the only possible exceptions are some oscillations in 6 Scuti 
stars 1 1 7 1 , PG 1159 stars 1 1 8 1 1 1 9 1 , and more recently, some very long periods of Rigel 1 20 1 , but none 
of these have yet been confirmed. In this connection, the identification of the e-mechanism as the 
driving mechanism behind the lightcurve variations of LS IV-14°116 would offer the first evidence 
that the e-mechanism can indeed drive pulsations in stars. 



2 Numerical methods and input physics 

Numerical models of ID-post-merger structures were constructed by simulating the merger process 
as an accretion process of the material of the secondary on top of the primary (a He-WD model). 
Specifically most of our models are constructed with MESA stellar evolution code |21| following the 
recipe detailed in |22| for fast merger episodes in which the accretion rate is of lO'^Mo/yr. Also some 
variants of the "composite merger" (see Il22l ) were computed. It is worth noting that nucleosynthesis is 
included during the accretion process. The composition of the accreted material was, by mass fraction, 
i2c=6xl0-^ "*N= 1.27x1 0-2, i6o=l.64xl0-\ 2"Ne=1.88xlO-^ and 24Mg=3.68xlO-^ for the Z = 
0.02 sequences and of l2c=3.1xlO-^ l4N^2.5xlO-^ "'0=9.3x10-^ 20Ne=3.9xlO-'* and24Mg=8.0x 
10 for the Z = 0.004 sequences. Hydrogen present in both the accretor and the accreted material was 
neglected. Evolution during the helium shell sub-flashes that follow the accretion process was followed 
with LPCODE stellar evolution code 1 23 1, which has been extensively used in the computation of He- 
flashes (e.g. 1241 and 125|). The pulsational instability analysis presented in this work was carried 
out with the hnear, nonradial, nonadiabatic pulsation code described in detail in |19| and references 
therein, in which the e-mechanism for driving pulsations in the He-burning regions is fully taken into 
account. As described in [15 | the "frozen-in convection" approximation was asssumed because the 
timescales of the pulsations are usually much shorter than the convective turnover times. 



2.1 Computed sequences 

In order to study the excitation of g-modes through the e-mechanism acting on the He-subflashes, and 
also to allow for comparisons with LS IV-14°116, we analysed the stability of modes with periods 
between 500s and 6000s. Modes were analysed every 5 stellar evolution timesteps during the whole 
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Table 1. Global properties of the merger sequences studied in this work compared with those of LS IV-14°116 
and the Hot-Flasher sequence of 1151 



Stellar Mass Initial Metallicity Range of H and He 





(Mo) 


(Zo) 


excited periods 


abundances by mass fraction 


0.25Mo+0.25Mq 


O.5OM0 


0.02 


< 3700s 


without H 


O.2OM0+O.2OM0 


O.4OM0 


0.02 


< 2600s 


without H 


0.25Mo+0.15Mo 


O.4OM0 


0.02 


< 2700s 


without H 


O.3OM0+O.IOM0 


O.4OM0 


0.02 


< 2500s 


without H 


0.25Mo+0.19Mo 


0.44Mo 


0.004 


< 2250s 


without H 


0.35Mo+0.09Mo 


0.44Mo 


0.004 


< 3440s 


without H 


0.25Mo+0.25Mo 


O.5OM0 


0.02 


< 3700s 


H=0.456 ,He=0.509 


0.25Mo+0.15Mo 


O.4OM0 


0.02 


< 2700s 


H=0.456 ,He=0.517 


Shallow Mixing fTSl 


0.47378Mo 


0.02 


< 2200s 


H=0.382 ,He=0.598 



LSIV-14°116[T3] - - 1954s-5084s H=0.544 ,He=0.452 



pre-EHB phase for all evolutionary sequences. The analysed sequences are detailed in Table[T]and the 
location of the sequences in the log g-\og Teff diagram is shown in Fig. [1] In addition, we analysed 
two sequences for which we added artificially a H-rich envelope with a H/He mixture similar to that 
inferred in LS IV-14°116. In order to maximize the effect of the presence of H, H was added as deep 
as possible without reaching layers with temperatures high enough to burn H (log(l - m{r)/Mstm) ~ 
-3.85). As will be shown, when comparing the log gf-log Teff of theoretical sequences and real stars, 
the total amount of H in the envelope should be taken into account. 



3 Discussion and Conclusions 

From Fig. [T]it is clear that more massive models evolve during the He-core flashes at higher tem- 
peratures than less massive ones. This is not unexpected as pure helium structures will settle on the 
so-called He main sequence ll26l which lies at higher temperatures for more massive models. Conse- 
quently the location of our sequences of higher masses is unable to fit the observed parameters of 
LS IV-14°116. It should be noted, however, that when H is added in the envelope of the theoretical 
sequences, tracks are shifted to lower temperatures and gravities, closer to the observed location of 
LS IV- 14° 116 (see top panels of Fig.[T]i. In the second and third row panels of Fig.[T] the evolution of 
the lower mass mergers analysed in this work is shown (Mflnai - O.4M0). As can be seen, lower mass 
models evolve in the pre-HB at lower temperatures (and through wider loops) than the more massive 
ones attaining temperatures similar to those of LS IV-14°1 16. In addition, it can be seen that the sur- 
face gravity of the sequences depend on how the final mass is attained, with mergers with initially 
more massive primary stars (i.e. lower accreted masses) evolving at lower gravities. 

As shown in Table[T]and Fig. |2]the e-mechanism in some of the more massive post-merger models 
is able to excite periods as long as 3700s, significantly longer than those predicted by the hot-flasher 
model studied in [ 15 |, making the predicted periods of the unstable modes closer to those of LS IV- 
14° 116. It should be noted however, that the formally unstable phase of the more massive sequence 
during the first flash (top panel. Fig. |2]i is too short, with the shortest e-folding times very similar to 
the length of the instability phase (~ lOOyr). Then, although that sequence displays formally unstable 
modes as long as 3700s they will not be observable. The situation is different for the less massive 
sequences (middle and bottom panels of Fig. |2]i, for which the instability phase is longer than the 
e-folding times, and thus pulsations driven by the e-mechanism should be observed. It is particularly 
worth noting that for the O.25M0+O.I5M0 sequence with H in the envelope (third row, right panel 
in Fig. [TJ, the remnant evolves through the observed location of LS IV-14°1 16 during the first three 
flashes, showing during this stage excited periods as long as the two shortest periods observed in LS 
IV-14°116 (P = 1953.74, 2620.27s|12|). In particular, the range of excited periods in the first five 
flashes is able to reproduce the higher amplitude mode of LS IV-14°116 (P - 1953.74). Still, the 
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Fig. 1. Location in the log g-log T^g of the post-merger sequences studied in this work. Thick color points indicate 
the location of models for which pulsations driven by the e-mechanism have been found with periods in the studied 
range. Color coding indicates the number of unstable modes found. Note that in all the sequences unstable modes 
are found during the fast redward evolution during the He-shell flashes. Cyan and grey boxes indicate the derived 
atmosphere parameters of LS IV-14°116 derived by L13J and L12J . respectively. 



excited modes of all our computed sequences have periods too short to explain the longest periods 
observed in LS IV-14°116 (P > 4000s). It should be noted, however, that the merger models we 
have adopted in this study are based on some simplifying assumptions, and thus more detailed merger 
models should be explored before making a final statement about the origin of LS IV-14°116 and 
the mechanism behind its pulsations. In particular, it is worth noting that the thermal structure of 
post-merger models is sensitive to the details of the merger process and on the relative masses of 
the primary and secondary stars [27 j . A more realistic treatment of the merger process might lead to 
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Fig. 2. Left: Evolution of tiie periods of the normal modes during the pre-EHB evolution. Excited periods through 
the e-mechanism are marlced with green thiclc points and compared with the observed modes of LS 1V-14°116 
(blue horizontal lines). Right: Same as left Fig.s but only during the first (and most violent) flash. Color coding 
indicates the value of the e-folding time of the unstable periods. The sequences plotted are the 0.25Mq + 0.25Mq 
merger (top), the 0.20Mo + 0.20Mo merger (middle) and the 0.25Mo + O.ISMq merger with an added hydrogen 
envelope (bottom). 



thermal structures signifficantly different and thus to a shift in the location of the He-core flashes, 
altering the range of excited periods. One interesting feature predicted by the present explanation of 
the pulsations observed in LS IV-14°116 is the large magnitude of the rate of period changes of the 
modes. In fact, as the structure changes very fast during the helium core flashes, then the predicted 
period change can be as high as a few seconds per year. The determination of large value in the rate of 
period changes would be a strong indication that LS 1V-14°1 16 is in a fast evolutionary stage. In this 
connection it might be worth noting that the value of the determined period of the ~ 2870s mode of 
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LS IV-14°116 differs in about 2s from the observations done in 2004 IfTOl and the determinations of 
2010 fT2\, which is somewhat larger than the quoted error bars (cr < Is). 

Despite the particular origin of LS IV- 14° 116, and it is by no means obvious whether it is a post- 
merger object or not, our analysis suggest that post HeWD-i-HeWD merger remnants should go through 
a phase in which pulations are driven by the e-mechanism. Due to the short duration of the unstable 
phases, only ~ 1 every 100 pre-EHB objects should be pulsators. A more detailed exploration of the 
parameter space of both the merger and hot- flasher scenarios to understand the pulsations excited in LS 
IV- 14° 1 16 will be performed in forthcoming works. We expect that future work will help to distinguish 
between the possible evolutionary channels for LS IV-14°1 16. 
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